First-principles study of high conductance DNA sequencing with carbon nanotube 

electrodes 
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Rapid and cost-effective DNA sequencing at the single nucleotide level might be achieved by mea- 
suring a transverse electronic current as single-stranded DNA is pulled through a nano-sized pore. 
In order to enhance the electronic coupling between the nucleotides and the electrodes and hence 
the current signals, we employ a pair of single-walled close-ended (6,6) carbon nanotubes (CNTs) as 
electrodes. We then investigate the electron transport properties of nucleotides sandwiched between 
such electrodes by using first-principles quantum transport theory. In particular we consider the 
extreme case where the separation between the electrodes is the smallest possible that still allows 
the DNA translocation. The benzene-like ring at the end cap of the CNT can strongly couple with 
the nucleobases and therefore both reduce conformational fluctuations and significantly improve the 
conductance. As such, when the electrodes are closely spaced, the nucleobases can pass through only 
with their base plane parallel to the plane of CNT end caps. The optimal molecular configurations, 
at which the nucleotides strongly couple to the CNTs, and which yield the largest transmission, 
are first identified. These correspond approximately to the lowest energy configurations. Then the 
electronic structures and the electron transport of these optimal configurations are analyzed. The 
typical tunneling currents are of the order of 50 nA for voltages up to 1 V. At higher bias, where 
resonant transport through the molecular states is possible, the current is of the order of several 
/j,A. Below 1 V the currents associated to the different nucleotides are consistently distinguishable, 
with adenine having the largest current, guanine the second-largest, cytosine the third and finally 
thymine the smallest. We further calculate the transmission coefficient profiles as the nucleotides 
are dragged along the DNA translocation path and investigate the effects of configurational varia- 
tions. Based on these results we propose a DNA sequencing protocol combining three possible data 
analysis strategies. 

PACS numbers: 87.14.G-, 87.80.St, 73.63.-b,73.40.Gk 



I. INTRODUCTION 

In 1996 Kasianowicz et al. demonstrated the possibility 
of translocating by electrophoresis single-stranded DNA 
(ssDNA) and RNA molecules through nano-sized biolog- 
ical channels^. This discovery initiated a surge of propos- 
als for novel DNA sequencing protocols 2 ^. Among the 
many Zwolak etal. proposed to sort the nucleobases from 
each other in a sequence by measuring the transverse tun- 
neling current across a ssDNA as this is driven through 
a nanoporei It was suggested that single base resolu- 
tion is in principle achievable, as the differences in the 
electronic structure among the nucleobases may translate 
into distinguishable tunneling current signals. However, 
as the tunneling process is extremely sensitive to the dis- 
tance between the electrodes and the nucleotides, one is 
left wondering whether the fluctuations in the transverse 
tunnel current, due to the fluctuations in the local geom- 
etry, may be too large to allow the identification of the 
various nucleotides. Furthermore additional fluctuations 
may arise from the electrostatic gate action introduced 
by the solvent molecules^. In order to overcome these 
intrinsic difficulties it was further suggested to trace the 
current fluctuations as the ssDNA is translocated across 
the pore, and to use their statistical distribution in order 
to unambiguously recognize the electronic signatures of 
the various nucleobases^—. 



Experimentally, Tsutsui et al. showed that it is pos- 
sible to identify single nucleotides in solution by two- 
probe tunneling current measurements and a thorough 
statistical analysis of the time-resolved current^. In this 
case the electrodes could be mechanically moved to op- 
timize the current through a given nucleotide. As such, 
although the experiment proves the concept of using elec- 
trical currents for the sequencing, it does not demon- 
strate a working device. 

Prototype sequencing devices with electrodes inte- 
grated in solid-state nanopores and nanofluidic channels 
have been also proposed^—. These have been realized 
experimentally by high resolution milling-based methods 
for a number of metal electrodes, but it has been specu- 
lated that similar techniques could be used for CNTs with 
possibly higher resolution^. On the theoretical front, 
the transverse tunneling conductance across nucleobases 
placed between two gold electrodes has been actively in- 
vestigated and debatec& 7 -'i£~— . Interestingly recently 
some special attention has been dedicated to exploring 
graphene nanoporc, graphene nanoribbon and carbon 
nanotubes (CNT) as potential electrodes materials^— . 

Despite these many works a key question still re- 
mains largely unanswered, namely how can one enhance 
the nucleotide-electrode interaction to a point where the 
transmigration is still possible, but the geometrical fluc- 
tuations are sufficiently suppressed to allow unambigu- 
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ous single nucleotide recognition. Possible strategics for 
achieving this goal are based on functionalizhig the elec- 
trodes with various chemical agents (including nucle- 
obases themselves for example) , which interact differently 
with the different nucclobase o 20 ' 24 ' 26 ' 27 . Yet, the func- 
tionalization may be extremely challenging in particular 
if this needs to be selective at the length scale of the typi- 
cal nanopore (a few nm). Here we use a different strategy 
that does not involve any chemical functionalization and 
considers single-walled close-ended (6,6) CNT electrodes 
placed at an ultra-short distance [see Fig. [T](a)]. An im- 
portant feature of these electrodes is that the benzene- 
like six-membered ring at the closed end cap provides 
the possibility of ir-ir coupling with the nucleobases, as 
ssDNA is transmigrated across. Importantly the ir-ir 
channel forces the nucleotides to align flat with respect 
to the electrode cap so that conformal fluctuations are 
highly suppressed. Our work analyses in detail the elec- 
tron transport in this favourable condition and propose 
a protocol for distinguishing the nucleotides. 

In short our computational strategy unfolds as follows. 
We first overview the general characteristics of the trans- 
mission coefficient as a function of energy for the four nu- 
cleotides sandwiched between two such electrodes. Then 
we search for the optimal molecular configurations of the 
nucleotides about the electrodes. For these configura- 
tions the electronic coupling is maximized and so is the 
low bias current. The electronic structure, the zero-bias 
transmission coefficients, and the I-V curves for all the 
nucleotides are analyzed next at the resulting optimal 
configurations. Finally we calculate the transmission co- 
efficient profiles along the full DNA translocation path 
and propose a DNA sequencing protocol combining mul- 
tiple data analysis. A method and a conclusion section 
complete this work. 



II. METHOD 

The device proposed here consists of a pair of semi- 
infinite close-ended CNT (6,6) electrodes aligned along 
the z-axis, and embedded inside a SisN4 nanopore [see 
Fig. [Q. The separation between the CNTs' caps is ini- 
tially only 6.6 A (this will be then optimized as described 
later). Since the carbon van der Waals radius is about 
1.7 A, one may expect that the small space between the 
electrodes will only accommodate the planar nucleobases 
with their base plane parallel to the plane of CNT end 
caps, but not the non-planar sugar ring and phosphate 
groups. This is indeed confirmed by molecular dynamics 
(MD) simulations performed with empirical potentials for 
a ssDNA pulled through a SiaN4 nanopore [see the snap- 
shot in Fig. Hlb)], from which we find that in general 
the nucleobases pass the gap between the CNTs approx- 
imately as indicated Fig. [Ha). A detailed description of 
the MD simulations for the DNA translocation is outside 
the scope of the present paper and it will be published 
elsewhere. Here we focus only on the transport proper- 




FIG. 1: The device setup considered in this work. In (a) 
we present the CNT/pXp/CNT (X=A, C, G, and T) junc- 
tion used to calculate the electron transport. The left- and 
right-hand side (6,6) CNT electrodes are semi-infinite and 
periodic along the transport direction (z-axis). In the scat- 
tering region the pXp molecules are sandwiched between the 
CNTs. A few layers of carbon atoms of both the CNTs are 
included in the scattering region allowing the convergence of 
charge density at the boundaries with the semi-infinite elec- 
trodes. In (b) a typical MD snapshot of the entire device, 
where the CNT electrodes are embedded into a SiaN4 -based 
nanopore [sketched by the circle in panel (a)], and an 11-base 
oligo-deoxyriboguanosine translocates through the nanopore 
(grey, cutaway view). The single-stranded DNAs translocate 
through the nanopore along the s-axis. A guanine base can 
be seen passing the gap between CNTs. Water and electrolyte 
ions are not displayed for clarity. 



ties of the ideal case. We note that not only does such a 
setup minimize the structural fluctuations during translo- 
cation, but it also leads to an effective 7r-7r bonding and to 
a strong electronic coupling between the nucleotides and 
both the electrodes. We therefore expect a much larger 
current than that for systems where the nucleotides arc 
oriented perpendicular to the electrodes, i.e. for the setup 
considered in several other studie s 20 ' 23 ' 24 . 

The four nucleotides commonly found in DNA are 
represented by four neutral chemical species, namely, 
adenosine 3'-phosphate-5'-phosphatc (pAp), cytidine 3'- 
phosphate-5 '-phosphate (pCp), guanosinc 3'-phosphate 
5'-phosphate (pGp) and thymidine 3'-phosphate-5'- 
phosphate (pTp). The two phosphate groups at both 
sides of each nucleobase are included in order to approx- 
imate the nearest chemical environment. In the initial 
configurations, each of the pXp (X=A, C, G, and T) 
molecules is sandwiched between the two CNT end caps, 
with their nucleobase mass centers coinciding with the 
origin. This is defined as the middle point of the elec- 
trode gap and the bonds connecting the sugar ring and 
the nucleobases aligned along the y-axis. We assume that 
the nucleotides translocate along the a;-axis. The unit 
cell size is set so that the distance between molecules in 
neighboring cells is at least 20 A, along both the x and 
y direction (we use periodic boundary conditions in the 
direction orthogonal to that of the transport, namely in 
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the x-y plane). 

The computation of the transmission coefficients (T) 
and the current-voltage (I-V) characteristics is per- 
formed by using the ab initio electronic transport 
code SmeagoIj2&~— . Smeagol implements the non- 
equilibrium Green's function (NEGF) method over den- 
sity functional theory (DFT), by using the pscudopoten- 
tial code Siesta^ as its electronic structure platform. In 
Smeagol the device under investigation is partitioned 
into three regions: the left- and right-hand side semi- 
infinite current- voltage electrodes, and the scattering re- 
gion [see Fig. [Ha)]. In our calculations the scattering re- 
gion comprises 156 and 144 carbon atoms respectively of 
left- and right-hand side CNT electrodes, and one of the 
pXp molecules. The portion of the electrodes included 
explicitly in the scattering region is large enough to en- 
sure that the calculated charge density at the outermost 
layers converges to that of bulk the CNTs. 

Both the electronic structure and the transport prop- 
erties are calculated with the local density approxima- 
tion (LDA) of the DFT exchange and correlation func- 
tional. A double-^ plus polarization basis sets is used for 
C(2s2p), N(2s2p), 0(2s2p), P(3s3p) and H(ls) through- 
out. All calculations are carried out with an equivalent 
real space mesh cutoff of 300 Ry and an electronic tem- 
perature of 300 K. In Smeagol the charge density is eval- 
uated by separating the integral over the non-equilibrium 
Green's function into an equilibrium part, which is per- 
formed along a contour in the complex energy plane, and 
a non-equilibrium part, which is performed along the real 
energy axis 2 ^ For the first wc use 16 energy points on 
the complex semi-circle, 16 points along the line parallel 
to the real axis and 16 poles. The integral over real en- 
ergies necessary at finite bias is evaluated over a mesh, 
whose energy spacing is not larger than 1 me V 28 ' 29 . 



III. RESULTS AND DISCUSSION 

A. General characteristics of the transmission at 
zero bias 

Wc first examine the general transport properties of all 
the possible four CNT/pXp/CNT (X=A, C, G, and T) 
junctions by plotting, in Fig. [3J the zero-bias transmission 
coefficient as function of energy for the initial geometries, 
where the nucleobases' center of mass coincides with the 
middle point of the electrode gap. Electron tunneling 
between the CNT electrodes (6.6 A apart in these initial 
simulations) through vacuum (we denote this configura- 
tion as CNT//CNT) is extremely small (T - 1.3 x 10~ 12 ) 
around the Fermi level, Ep. Note that the linear response 
(zero-bias) conductance, G, is simply G = Go T(Ep), 
where Go = 2e 2 /h is the quantum conductance, e is the 
electronic charge and h the Planck constant, so that a 
tiny T(E-p) is equivalent to a small linear response con- 
ductance. 

The inclusion of the pXp molecules in the device re- 




io~ 12 k 

-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 

E-E p (eV) 

FIG. 2: Zero-bias transmission coefficient as a function of en- 
ergy for CNT//CNT and the four CNT/pXp/CNT (X=A, C, 
G, and T) junctions. The curves refer to the initial config- 
urations where the nucleobase' center of mass coincides with 
the middle point of the electrode gap. The CNT electrodes 
are 6.6 A apart. 

suits in an increase of the transmission coefficient by over 
eight orders of magnitude with typical values of T(E-p) 
of around 10 _4 -10 -3 . We therefore expect in general 
a drastic change in the measured current when a nu- 
cleotide passes between the CNT electrodes, regardless 
of the nucleotide type. As such our device is at least 
capable of distinguishing whether or not a molecule is 
between the two electrodes. Notably the magnitude of 
the transmission coefficient at Ep is found to be orders 
of magnitude larger than that found by using either func- 
tionalized gold (~ 10~ 6 )2£ or graphene (> 10~ 10 )^ elec- 
trodes (in both cases the calculations have been carried 
out with Smeagol using similar computational param- 
eters). Clearly our much larger tunneling current origi- 
nates from the short electrode-electrode separation and 
by the fact that the nucleobases lie flat between the elec- 
trodes. 

In the various T(E) profiles there are peaks, as the one 
clearly visible at 0.5 eV, common to all the junctions. All 
of these are very pronounced and extremely sharp and 
they are characteristic of the CNT//CNT device as well 
(see peaks at -2.61 eV, -1.64 eV, -1.48 eV and 1.95 eV, 
with all the energies taken with respect to the CNT Ep). 
These are all due to localized surface states of the close- 
ended CNTs and therefore are not representative of any 
feature of the nucleotide electronic structure. Such CNT 
surface states can be readily identified by comparing the 
density of states (DOS) of an infinite CNT with the one 
of the semi-infinite close-ended CNT, shown in Fig. |3l 

For all the CNT/pXp/CNT junctions investigated we 
can identify two distinct transport regimes (Fig. [5]): 
a tunneling regime for energies between —1.0 eV and 
1.5 eV, and a resonant transport regime through the 
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FIG. 3: Density of states (DOS) of a (6,6) CNT elec- 
trode and of both the CNT//CNT and the CNT/pAp/CNT 
junctions as function of energy for the initial configura- 
tions. CNT/pAp/CNT is chosen as a representative of the 
CNT/pXp/CNT junctions. 
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FIG. 4: Total energy as function of the distance between the 
CNT electrodes, dcNT-CNT, for the CNT/pXp/CNT junc- 
tions (X=A, C, G, and T). The energy zero is taken at the 
lowest energy calculated. 



nucleotides' molecular levels for energies outside this 
range. Around E-p, the transmission curves for all the 
nucleotides are very smooth and vary little over energy. 
Importantly, the relative order of magnitude of the trans- 
mission coefficients, namely T p a p > T v q v > T p c p > TpT P , 
is preserved for a rather large energy range around the 
Fermi level. This indicates that it is possible to distin- 
guish the nucleotides already in the low-bias tunneling 
regime, and that the recognition is relatively robust with 
respect to fluctuations in the position of Ep. Further- 
more the flatness of the transmission coefficient around 
Ep suggests that one can infer the current for voltages 
up to about 1 V from the simple relation / w T(Ep)V . 
Therefore, in the next section we will analyze the trans- 
port properties at the Fermi energy only, and search 
for the molecular configurations of the CNT/pXp/CNT 
junctions that yield the largest transmission at the Fermi 
energy. When these configurations are explored by the 
ssDNA during the translocation process they will provide 
the largest contribution to the current . 



B. Search for the molecular configurations yielding 
the maximum transmission 

For the device setup investigated here we have iden- 
tified four possible degrees of freedom, which determine 
the relative position of the pXp molecules with respect 
to the electrodes: (1) the rotation angle of pXp around 
the z-axis [up to 60° due to the (6,6) CNT symmetry], 

(2) a translation of the pXp molecules about the y- and 

(3) the x-axis and (4) a variation of the distance between 
the two CNT electrodes. Here we assume that the inter- 
nal structure of the pXp molecules does not deform and 



that the nucleobases remain parallel to the CNT end cap 
planes during the translocation. 

In section UlI Fl we will discuss in details the effects of 
these operations on the electron transport; here we first 
identify the configurations which maximize the transmis- 
sion at E-p. To this goal we start from the initial con- 
figurations of the CNT/pXp/CNT junctions, defined in 
Sec. HIl and perform all the possible independent rota- 
tions and translations. Interestingly and importantly, we 
find that the molecular configurations having the largest 
transmission coefficients are generally found to be close 
to those which minimize the system energy. This is due 
to the fact that when the 7r-7r interaction between the 
nucleotides and the CNTs is large the total energy is 
generally minimized and the transmission maximized. 

We first determine the optimal distance between the 
CNTs. The guiding principle for the optimization is that 
the transmission is expected to increase for decreasing 
CNT separation, but the separation needs to be large 
enough to allow the translocation of the nucleotides. In 
any hypothetical fully operational device the CNT elec- 
trodes' position is fixed as these are embedded in the 
nanopore, so that we have to find a separation which is 
optimal for all the four nucleotides. In particular the 
electrodes gap should be wide enough to allow also the 
non-planar methyl group on thymine to pass through. 
In order to determine the optimal distance between the 
CNT electrodes we then fix the molecule in the lateral 
position (in the x-y plane) that yields the largest trans- 
mission at a CNT separation of 6.6 A and then change 
the separation between the electrodes, always keeping the 
molecule in the middle of the gap. 

We examine first the junction total energy as function 
of the separation between the CNT caps (see Fig.@|. The 
smallest possible CNT electrodes gaps can be seen to be 
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TABLE I: Optimal molecular configurations, which lead to 
maximal transmission at Ep, for the four CNT/pXp/CNT 
junctions (X=A, C, G, and T), relatively to the initial config- 
urations. 9 is the right-hand rotation angle of the nucleobases 
about the z-axis. The distance between the CNTs for optimal 
configurations is chosen to be commonly 6.4 A. 







pAp 


pCp 


pGp 


pTp 


9 




14.32° 


12.28° 


14.28° 


11.08° 


Offset 


»(A) 


-0.21 


-0.14 


-0.09 


0.03 


Offset 


x(A, peak 1) 


-0.05 


-2.13 


-0.58 


-1.88 


Offset 


a;(A, peak 2) 




1.01 




1.19 



approximately 5.8 A, since there is a sharp increase of the 
energy for closer separations regardless of the nucleotide 
type. Considering that the optimal distance may be un- 
derestimated by DFT calculations in which the van der 
Waals interactions (which tends to be repulsive at short 
distances) are not included, and that we should use a 
common distance for all four CNT/pXp/CNT junctions, 
in all subsequent calculations we assume that the device 
has been set up with a distance of 6.4 A between the two 
CNT electrodes. For this electrodes distance the molecu- 
lar configurations of the four CNT/pXp/CNT junctions 
yielding the maximum transmission at Ep are summa- 
rized in Table U in terms of the geometrical offset relative 
to their initial configurations. We denote these configu- 
rations as the optimal configurations. 

An interesting result of our search for the optimal con- 
figurations is that pCp and pTp turn out to have two con- 
figurations where the transmission coefficient peaks (see 
figure [5] for a graphic representation of all the optimal 
configurations). In these the smaller pyrimidine bases of 
the cytosine and thymine are located at either side of the 
CNT electrodes along the x-axis. In contrast, pAp and 
pGp appear to have only one major peak, and the config- 
uration is such that the larger purine bases of adenine and 
guanine are located in the center of electrodes gap. This 
difference between purine and pyrimidine bases results in 
distinct profiles of the transmission coefficient along the 
ssDNA translocation direction (see later in Section llll El) . 
In the calculations that will follow we will consider always 
the configuration responsible for the first of the transmis- 
sion peaks for the CNT/pCp/CNT and CNT/pTp/CNT 
junctions. 

C. Electronic structure and zero-bias transmission 
at the optimal configurations 

In order to understand the details of the electron trans- 
fer across the CNT/pXp/CNT (X=A, C, G and T) junc- 
tions, we now analyze the relation between the trans- 
mission and the electronic structure of the devices at 
their optimal configurations. In Fig. [B]wc plot the zero- 
bias transmission coefficient and the DOS of the pXp 
molecules forming the CNT/pXp/CNT junctions. The 
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c. CNT/pCp/CNT:peak 1 d. CNT/pCp/CNT:pi:ak 2 




c. CNT/pTp/CNT:pcak I f. CNT/pTp/CNTipeuk 2 




FIG. 5: Optimal molecular configurations for the 
CNT/pXp/CNT junctions (X=A, C, G, and T). For 
pCp and pTp the two equivalent configurations are both 
shown. 



vertical dashed lines indicate the energies of the eigenval- 
ues of the isolated pXp nucleotides, i.e. they correspond 
to the DOS of the isolated molecules. From the figure it is 
clear that the eigenvalues of the isolated molecules align 
well with the peaks in the DOS of the corresponding junc- 
tion (note that a global shift is applied in order to align 
the highest occupied molecular orbital -HOMO- of the 
molecule in the gas phase to that in the junciton). This 
is indicative of the fact that the molecule-electrode in- 
teraction does not involve charging, so that the molecule 
spectrum is not distorted in any significant way. How- 
ever the DOS of the junctions reveals a considerable level 
broadening, suggesting that there is a substantial overlap 
between the molecular orbitals of the nucleotides and the 
extended states of the electrodes. 

This strong overlap is responsible for the high trans- 
mission and Fig. [5] reveals a one to one correspondence 
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FIG. 6: Zero-bias transmission coefficient as function of en- 
ergy and corresponding density of states (DOS) projected 
onto the pXp molecules for the four CNT/pXp/CNT junc- 
tions. In the panels showing the DOS we also present, as 
dashed vertical lines, the DOS of the pXp molecules in the 
gas phase. 
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FIG. 7: Local density of state (LDOS) isosurfaces of the 
HOMO and LUMO of the CNT/pXp/CNT junctions at zero- 
bias. The geometries corresponds to the optimal configura- 
tions. 



TABLE II: Energy position of the HOMO and the LUMO 
of the pXp molecules in their corresponding CNT/pXp/CNT 
junctions (X=A, C, G, and T). 





pAp 


pCp 


pGp 


pTp 


HO (eV) 


-1.86 


-1.65 


-1.40 


-1.98 


LU (eV) 


2.06 


1.81 


2.42 


1.78 



between the molecule DOS and the peaks in the trans- 
mission functions. In particular for all the four junctions 
the HOMO and the lowest unoccupied molecular orbital 
(LUMO), both located on pXp molecules, are far away 
from the electrodes' Fermi level. As such the tunnel- 
ing transmission region around E-p corresponds to the 
molecule HOMO-LUMO gap. A summary of the energy 
position of both the HOMO and the LUMO of the four 
junctions is presented in Table ILTl 

At the optimal configurations (see Fig. [5J the pXp 
molecules are strongly coupled to the CNT electrodes, 
through the optimal ir-ir coupling of the nucleobases with 
the CNT end caps. In order to visualize such a 7r-7r bond, 
we calculate the real space local density of states (LDOS) 
associated to the charge density of both the HOMO and 
the LUMO of the CNT/pXp/CNT junctions. For the 
HOMO (LUMO) we integrate the charge density over an 
energy window 45 meV (50 meV) wide around the asso- 
ciated peak in the DOS. The results are shown in Fig. 
where the 7T-7T bonds are clearly visible. The face-to-face 
7r-7r stacking of the six-member rings of the CNT end 
caps with the nucleobases not only makes the optimal 
configurations energetically favorable, but also enhances 



the electron tunneling. Since the 7r-7r interaction disap- 
pears when the nucleobases are outside the central region 
in between the electrodes, one may expect that the dis- 
tribution of transmission coefficients (and hence the cur- 
rents) should peak when the nucleotides are inside such 
a region and that the high transmission configurations 
should be rather stable. 

Importantly the 7r-7r interaction of the pXp (X=A, C, 
G, and T) molecules with the CNTs involves only the nu- 
cleobases. This means that the sugar ring and two near- 
est phosphate groups in pXp do not play any significant 
role either in the bonding or in the electron transport. A 
demonstration of this feature is provided in Fig.[8j where 
we present the DOS projected respectively onto the nu- 
clcbase, the sugar and the phosphate groups. From the 
figure it is clear that both the HOMO and the LUMO are 
almost entirely associated to the nuclebases, while both 
the sugar and the phosphate groups contribute only to 
molecular levels far way from the Fermi level. This is 
a rather important aspect of the electronic structure of 
the DNA/CNT system investigated here, namely the fact 
that the non-sequence specific features of the DNA spec- 
trum are away from the electrodes Fermi level and there- 
fore are not expected to contribute significantly to the 
transport at least at relatively moderate voltages. 

Finally we wish to investigate further which functional 
part of the nucleotides contributes the most to the 7r-7r 
bonding by looking at the DOS projected either on the 
amine (NH 2 ) or the carbonyl (C=0) group (see Fig. [5]). 
Note that pAp, pCp and pGp have one primary amine 
group, pCp and pGp have one carbonyl and finally pTp 
has two carbonyls. From Fig.|H]we note that for pAp, the 
primary amine group contributes the most to the HOMO, 
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FIG. 8: Density of states (DOS) projected onto the base, the 
sugar and the phosphates of the various pXp molecules in 
the CNT/pXp/CNT junctions, calculated at zero-bias at the 
optimal configurations. 



with a major peak in the DOS at around —1.86 eV 
and minor peaks at around —1.66 eV and —1.44 eV. 
These form because of the coupling to the CNT elec- 
trodes. The presence of such minor peaks helps to ex- 
plain why at zero-bias pAp has the largest transmission 
coefficient despite its major HOMO peak lies below those 
of pCp and pGp (Table UTj) . In contrast for pCp, the car- 
bonyl group contributes more than the primary amine 
one to the HOMO, while for pGp these two groups con- 
tribute similarly to their HOMO. The main peaks of 
the carbonyl groups appear to be similar for pGp and 
pTp (around —2 eV) but these are located at higher en- 
ergies for pCp (around —1.65 eV). The different posi- 
tioning and strength of the electronic states located on 
these functional groups are important, since they affect 
the tunneling signals for the different nuclcobases. Their 
characteristic signals may be augmented by functionaliz- 
ing the electrodes with suitable chemical agents in order 
to amplify the electronic contrast between the different 
nucelobasea22»2&i2I. 



D. I-V characteristics at the optimal configurations 

We now investigate the transport properties of the pXp 
molecules (X=A, C, G, and T) by assuming that these are 
ideally measured at the optimal configurations. Fig. [TOl 
shows the zero-bias transmission coefficient, T(E), plot- 
ted on a logarithmic scale. As noted previously, the peaks 
in transmission at around 0.5 eV are due to localized sur- 
face states belonging to the CNT close-ends. These are 
unaffected by the pXp configuration and they contribute 
little to the electron current as bias is applied (for a dis- 
cussion on the resonant transport properties of surface 
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FIG. 9: Density of states (DOS) projected onto the amine 
and carbonyl groups of the various pXp molecules in the 
CNT/pXp/CNT junctions, calculated at zero-bias at the op- 
timal configurations. 
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FIG. 10: Zero-bias transmission coefficient as a function of 
energy for CNT/pXp/CNT (X=A, C, G, and T) junctions at 
the optimal configurations. 



localized state see reference [32]). At the Fermi level 
the relative order of magnitude of the transmission co- 
efficient is now TpAp > TpGp > T pCp > T pTp , i.e. it is 
the same order that we have found for the initial con- 
figurations (see Fig. [2]). However, if one now looks at 
energies away from E-p a few changes are notable. For 
instance at — 0.77 eV (energies are always taken from E-p) 
the transmission coefficient of pGp surpasses that of pAp 
and the order becomes T p g p > T p t p > T p a p > T p c p . In 
contrast for positive energies T p t p and T p c P first surpass 
T p Gp at 0.5 eV and then T p a p at 1.07 eV. In the region 
[1.1,1.8] eV, the order is now T pTp > T pCp > T pAp > 
T p Gpj which roughly reflects the energy order of the var- 
ious LUMOs (Table |TJ). Such a rather sensitive energy 
dependence of the order of the transmission coefficients 
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FIG. 11: I-V characteristics for the CNT/pXp/CNT (X=A, 
C, G and T) junctions at the optimal configurations in the 
[-1, 1] Volt bias range. 
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C, G and T) junctions at the optimal configurations in the 
[-4, 4] Volt bias range. 



for the various molecules reflects the fact that the nu- 
cleotides are electronically rather similar. However, we 
will show that the I-V curves offer a reasonable voltage 
range where the molecules can be distinguished. 

The calculated I-V curves in the moderate bias range 
[-1, 1] Volt are presented in FigfTT] In this particular bias 
window the transport is well within the tunneling regime, 
since there are no peaks in the transmission attributable 
to any pXp molecular orbitals (see Fig. ITU)) . Importantly 
in this moderate bias region the order of the currents am- 
plitude is the same as that of the transmission coefficients 
at the Fermi level, namely I p a p > / p g p > ^ P Cp > ^pTp- 
This is persistent throughout the entire bias values in- 
vestigated, suggesting that a current measurement in 
this range should be able to distinguish between the nu- 
cleotides. Furthermore, the current curves in FigfTT] are 
well spaced when the bias is between 0.7 V and 1 V, a bias 
interval which therefore emerges as the ideal probing volt- 
age range. We note that recently Tsutsui et al. reported 
a similar trend, namely Iq > Ic > It, for the peak 
currents of single nucleotides measured by two gold na- 
noclectrodes at a constant dc bias of 0.75 V—. 

In Fig. rT2] the I-V curves for the same configurations 
are shown over an extended bias window ranging now 
from —4 V to 4 V. From the figure we can clearly iden- 
tify the transition from tunneling to a resonant transport 
regime at about ±3 V for all nucleotides. This is the volt- 
age which coincides with a drastic increase of the current. 
For voltages larger than 3.5 V and smaller than —3.5 V, 
the current signals associated to the different nucleotides 
have a steady order, namely / pGp > / p t p > / P Cp ~ I P A P 
(Fig. [T2]) . Also notable is the fact that the I-V char- 
acteristics are approximately symmetric with respect to 
the bias polarity, since the nucleotides are placed in the 
middle of the gap between the CNT electrodes. 



E. Transmission profile of nucleotides along the 
translocation path 

In the discussion presented in the previous section, 
each of the four nucleotides is placed at the optimal con- 
figuration in the corresponding CNT/pXp/CNT (X=A, 
C, G, and T) junction. We now examine the evolution 
of the transport properties of the nucleotides when they 
translocate through the gap between the CNT electrodes 
along the x-axis. This is precisely the motion that the nu- 
cleotides will undergo in a nanopore experiment (Fig. [lj . 
For simplicity, the pXp molecules, which have limited 
conformational flexibility near the 6.4 A-wide electrode 
gap, are assumed to be rigid. Their positions are mea- 
sured as the distance between the nucleotide center of 
mass and the mid-point of CNT electrodes. 

Fig.fTSlshows the profiles of the transmission coefficient 
at the Fermi level, T(E-p) = Tq, along the nucleotides' 
translocation paths, with the peaks appearing at the op- 
timal configurations [see Table H] . As already mentioned 
in Section MI B| one can clearly see that the profiles of 
pCp and pTp are distinguishable from those of pAp and 
pGp by having one more pronounced transmission peak. 
The double peak structure arises from the two nearly 
degenerate optimal configurations, that both pAp and 
pGp possess because of their different symmetry. Since 
both the cytosine and the thymine base can be flipped 
before entering the electrode gap, the order and the rel- 
ative magnitude of these two peaks can be reversed. In 
addition we also note that the profile of pGp appears to 
have a broad shoulder, indicative of a second high trans- 
mission configuration along the translocation path. From 
Fig. [13] we also observe that the transmission coefficient 
decays sharply once the nucleotides is moved outside the 
electrodes gap. Notably the radius of the (6,6) CNT 
is about 4.1 A while the distance between adjacent nu- 
cleotides in a stretched ssDNA is about 7 A. As the tails 
of the transmission coefficient curves are quite small when 
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FIG. 13: Transmission coefficient at the Fermi level along the 
DNA translocation path. 

x < —3.5 A or x > 3.5 A (i.e., when the nucleobase is 
outside of electrodes gap) , we conclude that having nu- 
cleotides in a sequence should not affect much the reading 
of the tunneling current of a single nucleotide. 

We next integrate To along the translocation path to 
obtain the aggregate transmission coefficient (denoted as 
AT). Our results are summarized in Table Hill AT over 
the path -3.5 < x < 3.5 A accounts for 99 %, 94 %, 92 % 
and 96 % of the AT calculated over the more extended 
path —8 < x < 8 A, respectively for pAp, pCp, pGp 
and pTp. This is a further indication that most of the 
current is collected when the nucleobase is within the 
CNT region and that adjacent nuclcobases are far enough 
not to affect the measurement. Similar conclusions have 
been reached^ for Au electrodes at separations smaller 
than 1 nm. 

The aggregate transmission coefficient represents the 
time-average low-bias conductance over the translocation 
of each nucleotide, once we assume that the transloca- 
tion occurs at a constant velocity. Interestingly, we find 
a well spaced order in the AT, namely AT p a p > AT p q p > 
AT p c p > ATpT P , both when the AT is calculated within 
±3.5 A or ±8 A from the electrode gap center. Notably 
the order of the integrated transmission corresponds to 
the order found for the transmission at the optimal po- 
sitions, which indeed dominate the transport during the 
translocation. More importantly this fact indicates that 
the pAp>pGp>pCp>pTp order is recurrent in most of 
the transport measurements we have described and there- 
fore it is a rather robust result of the geometry investi- 
gated here. 



F. Effect of molecular configurations on the 
transport properties 

In the previous discussion about the transmission pro- 
files, we have assumed that the nucleotides translocate 



TABLE III: Aggregate transmission coefficient at the Fermi 
level along the DNA translocation path centered in the middle 
plane between the CNT electrodes (unit: 10 -4 A). 
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FIG. 14: Transmission coefficient at Fermi level for the 
CNT/pXp/CNT (X=A, C, G, and T) junctions as the nu- 
cleotides move along the y-axis. 

ideally along the x direction without any lateral displace- 
ment within the electrodes gap (along the y-axis). We 
now examine how Tq varies when the nuclcobases position 
fluctuates along the y direction, by using the nucleobase 
center of mass of the relative optimal configuration as ori- 
gin, y = 0. As it can be seen in Fig.[14l T decays quickly 
when the nucleobases move downwards away from the 
electrodes (negative y). This is similar to what we have 
observed in Fig.[13]for a longitudinal displacement. Note 
that there is not much room for positive y displacement 
now as this will need the non-planar sugar ring and phos- 
phate groups to enter the narrow electrodes gap. As such 
we stop the calculations at y = 1.5 A, which corresponds 
to the largest displacement possible along that direction. 

Next we consider the effect of rotations of the pXp 
molecules about the z-axis. As shown in Fig. [5l at the 
optimal configurations, the six-member ring of the end 
cap in each of the CNT electrodes is parallel to the plane 
of the nuclcobases. Moreover, the cap ring tends to align 
with the pyrimidine rings of the purine bases (adenine 
and guanine) so that the corresponding hexagon edges 
are parallel to each other. In contrast the edges arc stag- 
gered with the pyrimidine rings of cytosine and thymine 
so that one of its vertex atoms faces the center of the 
pyrimidine. One may then expect that rotations may 
produce significant change in the low-bias conductance. 
This is however not the case. In fact by plotting Tq as a 
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FIG. 15: Transmission coefficient calculated at the Fermi level FIG. 16: Transmission coefficient at Fermi level for the 

for the CNT/pXp/CNT (X=A, C, G, and T) junctions as a CNT/pXp/CNT (X=A, C, G, and T) junctions as a func- 

function of the rotation of angle of the pXp molecules about tion of the electrode separation, dcNT-CNT- Note that To 

the z-axis (vertical). decays exponentially when dcNT-CNT increases. 



function of the rotation angle, 9, (Fig. [T5|) we find only 
small changes, and most importantly also that the order 
of the amplitudes of the conductance of the various nu- 
cleobases does not change. This essentially means that 
rotations about the optimal configurations will not un- 
dermine our ability to identify the nucleobases. 

Finally, we investigate the effects that the gap size be- 
tween the CNT electrodes has on the transport by plot- 
ting Tq as function of the electrode distance, c?cnt-cnt 
(Fig. n~6|) . In this case we maintain the position of the 
plane containing the pXp molecules fixed at the mid- 
point between the electrodes. As expected from a tun- 
neling process we find that the transmission coefficient 
decreases exponentially with increasing <icNT— cnt- The 
figure [16] also shows that the exponential decay rate a, 
defined as the slope of the lnT(d) curves, is rather similar 
among the nucleotides. Interestingly To for pAp is always 
well separated from that of the other nucleotides for all 
the electrode gap distances investigated. In contrast the 
relative order between pCp, pGp and pTp depends on 
^CNT-CNT- It then emerges that it is absolutely crucial 
to maintain the electrode separation as stable as possi- 
ble during a measurement and also that this should be 
the smallest possible, so that configurational fluctuations 
will be highly suppressed. 

G. DNA sequencing protocol 

With all the results of the previous sections in mind 
we can now formulate a possible protocol for reading the 
nucleotide sequence. The proposed device functions by 
translocating ssDNA through a nanopore and simultane- 
ously by measuring the transverse current between two 
closely spaced CNT nanotubes [(6,6) close-ended CNTs 
at 6.4 A separation in our case]. The DNA strand can 
be pulled through the nanopore mechanically by an opti- 



cal or magnetic tweeze r i with an approximately con- 
stant velocity. Typical velocities for optical tweezers have 
been reported^ at around 30 rrm-s . This means that 
one need a time resolution in the kHz range in order to 
obtain a spatial resolution of 0.1 A, which is certainly 
enough for detecting accurately the transmission profiles 
of the nucleotides shown in Fig. [T31 Note that the cur- 
rents involved are of the order of 10 nA and we need an 
amplitude resolution of the order of 1 nA or less (de- 
pending on the bias) for typical resistances of the order 
of 10 Mf2. These are rather standard demands for con- 
ventional electronics. 

In contrast if the ssDNA translocation is electrophoret- 
ically driven, then the velocity may fluctuate consider- 
ably so that the recognition of the current profiles may 
become difficult. Furthermore the typical translocation 
velocities are much larger than those achievable by me- 
chanical manipulations. For instance early studies^ re- 
ported speeds of the order of 10 6 nm-s _1 , which require 
GHz time-resolution for a sub-A spatial resolution. More 
recently a substantial slowing down has been achieved 
by gate modulation of the nanopore wall surface charges, 
with reported velocities of 55 ^m-s -1 . This translates 
in a 5 MHz time resolution for 0.1 A spatial resolution. 
Clearly GHz frequencies are out of the reach, unless the 
typical resistances are in the kfi range, while MHz reso- 
lution may be possible. In any case we believe that for 
a usable device setup electrophoretically driven translo- 
cation is not suitable at present, at least if the signal is 
made of the current profiles during translocation. 

Further information however can be extracted by sim- 
ply taking the integral over time of each current peak 
and/or successive peaks. These are acquired during the 
translocation of an entire nucleobase across the gap be- 
tween the electrodes. As such the time-resolution re- 
quired is now about two orders of magnitude larger than 
that needed for resolving the profiles. This means oper- 
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FIG. 17: Transmission coefficient at the Fermi level against 
the interaction energy for the CNT/pXp/CNT (X=A, C, G, 
and T) junctions along DNA translocation path. 



ating the device at ^100 kHz frequencies for slow elec- 
trophoretically driven translocations. 

Finally an alternative quantity to measure is, as pro- 
posed first by Zwolak et al.—, the statistical distribution 
of the current over time. This is a measure probably 
achievable only with the tweezers setup and essentially 
consists in holding the ssDNA at a fixed position and 
then, after the system has thermally equilibrated inside 
the nanopore, in measuring the current. We expect that 
the details of the current distributions may be device- 
dependent, but the peak currents of the nucleobases 
should maintain the order / p Ap > IpGp > ^pCp > IpTp, 
at least in the suitable low bias range around [0.7, 1] V. 

The statistical distributions of the transmission coeffi- 
cient/current of the CNT/pXp/CNT junctions are yet to 
be calculated for our device (Fig. [1]). However, some in- 
sights can be obtained by plotting the To together with 
the interaction energy between the pXp molecules and 
the CNT electrodes along the translocation path. The 
interaction energy is defined as AE = E p x P + Eqnt — 
-E'CNT/pXp/CNT: where E a is the DFT total energy of 
the a system (a = pXp, CNT and CNT/pXp/CNT). 
Our results are shown in Fig. [T7] The calculated inter- 
action energies are one order of magnitude larger than 
that for typical hydrogen bonding (0.1 — 0.4 eV). This 
is due to the strong ir-ir interaction found in our device 
at the optimal configurations. Interestingly, we find that 
the position of the maximum transmission corresponds 
approximately to the position displaying the largest in- 
teraction energy. As such, one may expect that when a 
nucleobase passes through the device, it will tend to re- 
main between the electrodes and parallel to the end caps 
of CNTs. In other words the same interaction responsible 
for the relatively large currents and for the suppression 
of the geometry fluctuations also acts to slow down the 
DNA translocation across the nanopore. 



IV. CONCLUSIONS 

In summary we have proposed a device for high 
throughput DNA nanopore sequencing, based on reading 
the transverse current across the nucleobases of ssDNA as 
it passes between two closely spaced single- walled close- 
ended (6,6) CNTs. The device operation has been inves- 
tigated theoretically with state of the art density func- 
tional theory combined with non-equilibrium quantum 
transport. In general the electron transport through the 
various nucleobases can be roughly separated into a tun- 
neling regime at low bias and a resonant transport one at 
higher bias. We have identified the optimal position of 
the various nucleotides yielding the highest zero-bias con- 
ductance and used them to analyze the various parame- 
ters affecting the device. In particular we have focussed 
our attention on the dependence of the electron transmis- 
sion on the longitudinal position of the nucleotides with 
respect to the center of the electrodes, on possible lateral 
displacements and on the electrode separation. 

The current/conductance profiles of the nucleotides 
along the DNA translocation path reveal that the pyrimi- 
dine bases (cytosine and thymine) have a two peaks struc- 
ture, indicative of two high transmission geometries. The 
same is not found for the other two bases. Importantly 
we have found that the peak tunneling currents of the nu- 
cleobases have the following order / p a p > / P Gp > ^pCp > 
JpTp at low bias and that the same order is maintained 
for the aggregate transmission coefficients. Furthermore 
we have demonstrated that at this close electrode sep- 
aration the current drops drastically as the nucleotides 
exit the electrodes region, so that well separated current 
signals should be detected during the translocation. 

Finally we have proposed a measurement protocol for 
rapid DNA sequencing. This is presently suitable for 
tweezer driven translocation and consists of three inde- 
pendent measurements. The first involves the pattern 
recognition of the current profile during the transloca- 
tion of the nucleobases. The second relies on the fact 
that the aggregate transmission coefficient (and hence 
the integral of the current) along the translocation path 
has a well-defined order. Finally the third concerns the 
time-distribution of the current for mechanically stabi- 
lized junctions, where the nucleotide to measure is held 
at a fixed position. 
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